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Abstract: 

We have calculated the complete electroweak 0(a) radiative corrections to the single 
Higgs-boson production processes e + e~ — > z/^H (I = e, /x, r) in the electroweak Standard 
Model. Initial-state radiation beyond 0(a) is included in the structure-function approach. 
The calculation of the corrections is briefly described, and numerical results are presented 
for the total cross section. In the scheme, the bulk of the corrections is due to initial- 
state radiation, which affects the cross section at the level of —7% at high energies and 
even more in the ZH threshold region. The remaining bosonic and fermionic corrections 
are at the level of a few per cent. The confusing situation in the literature regarding 
differing results for the fermionic corrections to this process is clarified. 
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1 Introduction 



The investigation of the mechanism of electroweak symmetry breaking in general and 
of the Higgs boson in particular will be one of the main tasks at future colliders. While the 
LHC will discover the Higgs boson, if it exists and has no particularly exotic properties, 
its complete profile can only be studied in the clean environment of an electron-positron 
linear collider. These studies require adequate theoretical predictions including radiative 
corrections and finite-width effects. 

In e + e~ annihilation there are two main production mechanisms for the Standard 
Model (SM) Higgs boson. The cross section of the Higgs-strahlung process, e + e~ — > ZH, 
rises sharply at threshold to a maximum a few tens of GeV above the threshold energy 
M z + M H and then falls off as where ^Js is the centre-of-mass (CM) energy of the 
e + e~ system. In the W-boson fusion process, e + e~ — > z/ e z/ e H, the incoming e + and e~ each 
emit a virtual W boson which fuse into a Higgs boson. The corresponding cross section 
grows as In s and thus is the dominant production mechanism for large energies. 

In lowest order, the Higgs-strahlung process has been studied in Ref. [Q] and the 
vector-boson fusion process in Ref. [E]. The O(o) electroweak radiative corrections to 
the process e + e~ — > ZH have been calculated by different groups [ Ej many years ago. 
The electroweak corrections to e + e~ — > vvH have attracted a lot of interest recently. 
In Ref. [0j the leading corrections in the limit of a heavy top quark, as well as the 
corresponding QCD corrections in 0(a^G^m1) with n = 1,2, have been worked out. The 
contributions of fermion and sfermion loops in the Minimal Supersymmetric Standard 
Model have been evaluated in Refs. [EH EH; at first sight, however, the results of the two 
calculations do not agree, not even on the fermion-loop contributions in the SM. A first 
calculation of the complete 0(a) electroweak corrections to e + e~ — > vuH in the SM has 
been performed very recently [|7j. Analytical results for the one-loop corrections to this 
process have also been obtained by another group [IHj as MAPLE output, but a numerical 
evaluation of these results is not yet available. 

In this paper we present first results of a completely independent calculation of the 
0(a) electroweak corrections to the complete process e + e~ — > vvH in the SM. Details on 
this calculation will be given elsewhere. Here we sketch only the main ingredients. 

2 Method of calculation 

We have calculated the complete 0(a) electroweak virtual and real photonic correc- 
tions to the processes e + e~ — > u e u e B., u^P^H, and v T v T H. For e + e~ — > z/ e z/ e H, this includes 
both the corrections to the Higgs-strahlung and the vector-boson fusion processes, which 
are taken into account coherently. 

The calculation of the one-loop diagrams has been performed in the 't Hooft-Feynman 
gauge both in the conventional and in the background-field formalism using the conven- 
tions of Refs. [E] and [HD1> respectively. The renormalization is carried out in the on-shell 
renormalization scheme, as described there. The electron mass m e is neglected whenever 
possible. 

The calculation of the Feynman diagrams has been performed in two completely inde- 
pendent ways, leading to two independent computer codes for the numerical evaluation. 
Both calculations are based on the methods described in Ref. [ E] . The tensor coefficients 
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of the one-loop integrals are recursively reduced to scalar integrals with the Passarino- 
Veltman algorithm [[IT] at the numerical level. The scalar integrals are evaluated using 
the methods and results of Refs. [E1H2], where ultraviolet divergences are regulated di- 
mensionally and IR divergences with an infinitesimal photon mass. The two calculations 
differ in the following points. In the first calculation, the Feynman graphs are generated 
with FeynArts version 1.0 [U^J. Using Mathematica the amplitudes are expressed in terms 
of standard matrix elements and coefficients of tensor integrals. Tensor 5-point functions 
have been evaluated both by applying the usual Passarino-Veltman reduction and by us- 
ing the direct reduction to 4-point integrals of Ref. [Ej. While the results based on the 
Passarino-Veltman algorithm become numerically unstable at the phase-space boundary 
owing to the appearance of inverse Gram determinants and could only be rescued by a 
careful extrapolation out of the numerically safe inner phase-space domains, the direct 
reduction of Ref. avoids inverse leading Gram determinants, rendering the results of 
this approach well behaved near the phase-space boundary. The whole calculation has 
been carried out in the conventional and in the background-field formalism. The second 
calculation has been done with the help of FeynArts version 3 [E3 and FormCalc [ ITS]. 
The analytical expressions generated by FormCalc were translated into C code. In or- 
der to eliminate 5-point tensor integrals the interference of the pentagon diagrams with 
the lowest-order amplitude was evaluated with FeynCalc [ IT7j by evaluating the fermion 
traces. Then the loop momenta of the 5-point integrals appeared in the numerator only in 
scalar products which could be cancelled with denominators, leaving only scalar 5-point 
integrals. 

The results of the two different codes, those obtained within the conventional and 
background-field formalism, and those resulting from different treatments of the tensor 
5-point functions are all in good numerical agreement (typically within at least 12 digits 
for non-exceptional phase-space points). 

We use two different schemes for the inclusion of the finite Z-boson decay width. In 
the fixed-width scheme, each resonant Z-boson propagator l/(s uP — M|), where s uP is the 
invariant mass of the neutrino-antineutrino pair, is replaced by l/(s vD — M| + iM^Ti), 
while non-resonant contributions are kept untouched. As a second option, we applied a 
factorization scheme where the full (gauge-invariant) ZH-production amplitude with zero 
Z-boson width is rescaled by a factor (s u p — Mf )/ (s ut > — M| + iMzPz) • Within integration 
errors both schemes give the same results for the total cross section. 

The matrix elements for the real photonic corrections are evaluated using the Weyl- 
van der Waerden spinor technique as formulated in Ref. [E] and have been success- 
fully checked against the result obtained with the package Madgraph [d]. The soft 
and collinear singularities are treated both in the dipole subtraction method following 
Refs. [I2H1I2J] and in the phase-space slicing method following closely Ref. [122]- Beyond 
0(a) initial-state-radiation (ISR) corrections are included at the leading-logarithmic level 
using the structure functions given in Ref. [|2H] (for the original papers see references 
therein) . 

The phase-space integration is performed with Monte Carlo techniques in both com- 
puter codes. The first code employs a multi-channel Monte Carlo generator similar to 
the one implemented in RacoonWW [IIHJ 121] and Lusifer [I2H], the second one uses the 
adaptive multi-dimensional integration program VEGAS [125]. 



2 



3 Numerical results 



3.1 Input parameters 

For the numerical evaluation we use the following set of SM parameters [I27j. 





= 1.16639 x l(r 5 GeV- 2 , 


or(0) 


= 1/137.03599976, 






M w 


= 80.423 GeV, 


M Z LEP 


= 91.1876 GeV, 


pLEP 
1 Z 


= 2.4952 GeV, 


m e 


= 0.510998902 MeV, 


"fy 


= 105.658357 MeV, 


m T 


= 1.77699 GeV, 


m u 


= 66 MeV, 


m c 


= 1.2 GeV, 


m t 


= 174.3 GeV, 


m d 


= 66 MeV, 


m s 


= 150 MeV, 


m h 


= 4.3 GeV. 



(3-1) 

We do not calculate the W-boson mass from G M but use its experimental value as input. 
Since we employ a fixed width in the resonant Z-boson propagator in contrast to the 
approach used at LEP to fit the Z resonance, where a running width is taken, we have to 
convert the "on-shell" values of M Z EP and T^ EP , resulting from LEP, to the "pole values" 
denoted by Mz and Tz in this paper. The relation of the two sets of values is given by [ 

EH! 

M z = M% EP /yJl + (r£ EP /Af| J3P ) 2 = 91.1535 GeV, 

r z = r£ EP / yfl + (rLEp/M^p) 2 = 2.4943 GeV, (3.2) 

i.e. the difference is of formal two-loop order and numerically hardly visible in the results 
presented below. The masses of the light quarks are adjusted to reproduce the hadronic 
contribution to the photonic vacuum polarization of Ref . [ 29] . Since we parametrize the 
lowest-order cross section with the Fermi constant G^ (G^ scheme), i.e. we derive the 
electromagnetic coupling a according to a^^ = \f2G tl My f s 2 v /'K, the results are practically 
independent of the masses of the light quarks. Moreover, this procedure absorbs the 
corrections proportional to m 2 /M^ in the fermion- W-boson couplings and the running 
of ct(Q 2 ) from Q 2 = to the electroweak scale. In the relative radiative corrections, we 
use, however, a(0) as coupling parameter, which is the correct effective coupling for real 
photon emission. 

We always sum over all three neutrino species, i.e. over the processes e + e~ — > z/ e z/ e H, 
v^jp^H, and u T u T H. Besides the full cross section, denoted "total" in the plots, we also 
give the cross section resulting from the ZH-production channel and the WW-fusion chan- 
nel separately, which are referred to as "ZH" and "WW" contributions, respectively. In 
the ZH-production channel we sum over the relevant contributions of all z/z/H final states, 
which is equivalent to multiplying the cross sections for e + e~ — > v^v^Yi by a factor 3. This 
means that the results shown for "total" and "ZH+WW" only differ by the interference 
terms between the ZH and WW channels. In the results presented here, the ISR is con- 
voluted only with the lowest-order cross section. We consider merely total cross sections 
without any cuts; distributions will be discussed elsewhere. 

For reference we give some numbers for the total cross section in lowest order, a tTee , 
and including electroweak corrections, a, together with the relative corrections defined as 
S = a/a tTee — 1 in Tabled] The last numbers in parentheses correspond to the Monte 
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M H [GeV] 


O-tree [fb] 


a [fb] 


S[%] 


115 


92.64(2) 


85.01(8) 


-8.2(1) 


150 


68.17(2) 


61.76(5) 


-9.4(1) 


200 


41.800(9) 


37.76(3) 


-9.7(1) 


250 


23.764(4) 


20.97(1) 


-11.8(1) 


300 


12.125(2) 


10.478(6) 


-13.6(1) 


350 


5.2047(6) 


4.264(2) 


-18.1(1) 



Table 1: Lowest-order cross section for e 4 
tion including full electroweak corrections, 
masses at y/s = 500 GeV 

Otree [fb] 
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Figure 1: Lowest-order cross section and contributions from ZH-production and WW- 
fusion channels for Mh = 115 GeV and Mh = 150 GeV 



Carlo integration error of the last given digit. 

In Figure Q we show the lowest-order cross section as a function of the CM energy for 
Mh = 115 GeV and 150 GeV. Besides the total contribution we separately give also the 
contributions from ZH production and WW fusion as defined above. Below the ZH thresh- 
old and above 400 GeV the cross section is dominated by WW fusion. The ZH-production 
contribution dominates from the ZH threshold up to about 300 GeV. Comparing the "to- 
tal" cross section with the sum "ZH+WW" , one can see that the interference contributions 
between ZH and WW channels are small in lowest order for the inspected kinematical 
situation. 

The relative corrections to the lowest-order contributions of Figure ^ are shown in 
Figure 121 The corrections to the ZH-production channel are large and negative (< —20%) 
below threshold, rise fast above threshold, and reach 19% and 11% at y/s = 1 TeV for 
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Figure 2: Relative electroweak corrections to the complete process e + e — ■> vvH and to 
the contributions from ZH-production and WW-fusion channels for Mh = 115 GeV and 
M H = 150 GeV 



M H = 115 GeV and M H = 150 GeV, respectively. The corrections to the WW-fusion 
channel are similar below the ZH threshold, rise sharply at the threshold and are flat and 
about —10% above 500 GeV. The corrections to the complete process are always negative. 
They follow those of the ZH-production channel for energies below ~ 300 GeV and those 
of the WW-fusion channel above ~ 500 GeV. 

The contributions of 0(a) ISR corrections, ISR beyond 0(a), fermionic corrections, 
and non-ISR bosonic corrections to the ZH-production and WW-fusion channels are pre- 
sented in Figure El The 0(a) ISR corrections are defined as the O(o) contributions of 
the structure functions of Ref. [E3j- The fermionic corrections summarize all fermion-loop 
contributions to loop diagrams and counter terms. The non-ISR bosonic corrections are 
obtained by subtracting the 0(a) ISR corrections and the fermionic corrections from the 
complete O(a) corrections. It can be seen, that the energy dependence in both channels 
and, in particular, the rise of the corrections to the ZH-production channel with energy 
is essentially due to the ISR corrections. The large corrections above the peak of the 
lowest-order ZH cross section are due to the decreasing cross section. Since the ISR cor- 
rections increase with the height of this peak, they are larger for smaller Higgs masses (see 
Figure I2J). The ISR corrections beyond 0(a) are at the level of several per cent where the 
lowest-order cross section rises strongly, but are below 1% above 300 GeV. The non-ISR 
corrections provide a measure of the genuinely weak corrections. For the ZH-production 
channel (Figure El left) the fermionic corrections increase slowly from 3% to 11% around 
500 GeV and then go down to 9% at 1 TeV. The non-ISR bosonic corrections decrease 
from —5% to —20% over the considered energy range. This behaviour is typical if elec- 
troweak Sudakov logarithms of the form —a In 2 (s/M|) dominate the corrections, which 
is always the case if the major part of the total cross section results from intermediate 
scattering angles. For the WW-fusion channel (Figure E3 right), the ISR corrections are 
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Figure 3: Relative electroweak corrections to the ZH-production (left) and WW-fusion 
(right) channels resulting from 0(a) ISR, ISR beyond 0(a), fermion loops, and non-ISR 
bosonic corrections for Mh = 150 GeV 
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Figure 4: Relative electroweak corrections to the complete lowest-order cross sections re- 
sulting from 0(ot) ISR, ISR beyond 0(a), fermion loops, and non-ISR bosonic corrections 
for M H = 150 GeV 
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negative for all energies, which can be explained by the rising cross section with increasing 
scattering energy. ISR beyond 0(a) influences the cross section only at the level of < 1% 
for energies above 300 GeV. The fermionic corrections are between —1% and —2% over 
the whole energy range. The non-ISR bosonic corrections fall from +1% at threshold to 
—3% at a/s = 1 TeV. This energy dependence is much weaker than for the ZH chan- 
nel, since the cross section for the WW channel is more and more dominated by small 
scattering angles, i.e. not by the Sudakov regime. 

The various contributions of the corrections to the complete process are depicted in 
Figure 0] The size of the corrections can be easily explained by the results of Figure H2 
it is given by the size of the corrections of the dominating channel. The ISR corrections 
vary strongly in the region of the ZH threshold but are nearly flat and about —7% for 
energies above 400 GeV. They are always negative since the lowest-order cross section is 
continuously rising. The fermionic corrections reach a maximum of about 6% in the region 
where the ZH-production channel dominates and are between 0% and — 1% above 500 GeV 
where the WW-fusion process dominates. The non-ISR bosonic corrections exhibit a 
minimum of about —7% where the ZH-production channel dominates and are between 
0% and —4% elsewhere. 

4 Comparison with other calculations 

We have compared our results for the 0(a) corrections to Ref. [Ej and the contribu- 
tions from closed fermion loops with Refs. [EIIE1- 

Adapting the input parameters and the parametrization of the lowest-order matrix 
element to those used in Ref. [|7j, we reproduced the numbers in Table 2 for the total 
cross section given in the first paper of Ref. [|7j. We find agreement within 0.2% for 
the total lowest-order cross section 1 and within 0.3% for the corrected cross section. 
The corrections relative to the lowest-order cross section agree within 0.2%. Note that 
the authors of Ref. [Ej use a(0) to parametrize the lowest-order cross section. As a 
consequence their relative corrections are shifted by 3Ar « +9% compared to those in 
the scheme. 

Adapting the input parameters to those of Ref. [E], who also use the scheme, we 
perfectly reproduce the SM results of these authors for the tree-level cross section and the 
fermionic corrections. 

The authors of Ref. use a(Mz), i.e. the running electromagnetic coupling in the 
MS scheme at M z , as defined in Eq. (B.2) of Ref. [EU], to parametrize the lowest-order 
cross section. The relative fermionic corrections in this scheme differ from those in the G M 
scheme by 3[(Ar) ferm — Ad(M z )] « —13%. Moreover, Eberl et al. take into account only 
the loops from the top-bottom doublet and omit the corrections to the ZH-production 
channel. Adapting to this setup, we find agreement with the results of Ref. at 1 TeV 
within the integration errors for the lowest-order cross section. The corrections relative 
to the lowest-order cross section agree typically within 0.3%. The large differences in 
relative corrections between Refs. [ |5] and [ E] thus result essentially from the use of 

*As we were told by the F. Boudjema, the integration errors of the numbers for the lowest-order cross 
section in Table 2 of Ref. [[7], which were suppressed in the table, are also of the order of 0.2%. When 
increasing statistics the agreement becomes better than 10~ 4 for the lowest-order cross section. 
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different parametrizations of the lowest-order matrix element. The fact that the lowest- 
order SM cross sections in Refs. and [E] agree qualitatively despite of the different 
input-parameter schemes used is accidental and due to the (different) input parameters. 

Where the WW-fusion channel dominates, the non-ISR, i.e. the genuine weak, correc- 
tions are large if a(0) or cx(Mz) are used to parametrize the lowest-order cross section. 
Only when using the large corrections associated with the running of a and those pro- 
portional to m\ in the W-boson-fermion couplings are absorbed in the lowest-order cross 
section. For the ZH-production channel the situation is more complicated and a dedicated 
improved Born approximation is under investigation. It should be noted that higher-loop 
corrections become relevant in schemes where the one-loop corrections are large. Since 
ISR corrections, the running of a, and the m t -dependent corrections are independent, one 
has to deal with all these terms separately. 

5 Summary 

We have presented results from a calculation of the complete electroweak 0{a) ra- 
diative corrections to the single Higgs-boson production process e + e~ — > uuB. in the 
electroweak Standard Model. We find that the ISR corrections are of the order of —7% 
at high energies and more than —10% near the ZH threshold. In this region even the 
higher-order ISR corrections reach several per cent. This is due to the strong energy 
dependence of the lowest-order cross section. The non-ISR corrections are at the level of 
a few per cent if the lowest order matrix element is parametrized with the Fermi constant 
Gfj,. In other schemes these corrections are of the order of 10%. It has been pointed out 
that the confusion in the literature regarding the size of the corrections to e + e~ — > vvH 
is due to the use of different schemes and input parameters. 
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